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Abstract I summarize recent results, obtained with E. Demler, K. Park, A. Pol- 
kovnikov, M. Vojta, and Y. Zhang, on spin and charge correlations near 
a magnetic quantum phase transition in the cuprates. STM experiments 
on slightly overdoped Bi2Sr2CaCu208+i (J. E. Hoffman et al., Science 
295, 466 (2002)) are consistent with the nucleation of static charge 
order coexisting with dynamic spin correlations around vortices, and 
neutron scattering experiments have measured the magnetic field depen- 
dence of static spin order in the underdoped regime in La2_5SraCu04 
(B. Lake et al, Nature 415, 299 (2002)) and LaCu04+y (B. Khaykovich 
et al, Phys. Rev. B 66, 014528 (2002) ). Our predictions provide a 
semi-quantitative description of these observations, with only a single 
parameter measuring distance from the quantum critical point chang- 
ing with doping level. These results suggest that a common theory of 
competing spin, charge and superconducting orders provides a unified 
description of all the cuprates. 
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Introduction 

Three recent experiments [1, 2, 3] have shed new hght on the spin 
and charge density wave collective modes of the cuprate superconduc- 
tors. This article will summarize the main results of our theory [4, 5, 
6, 7, 8, 9, 10] of these collective modes in the vicinity of a quantum 
phase transition between two superconducting states, only one of which 
has static, long-range, spin density wave order: we will recall our main 
predictions, and discuss further experimental tests, wc will also connect 
our theory to these experiments. In particular, we suggested [6] that 
static charge order should coexist with dynamic spin-gap fluctuations 
around the vortex cores in the cuprate superconductors, as may have 
been been observed in [1]. Wc also discussed [5] a singular field depen- 
dence for the static magnetic moment in the underdoped cuprates, and 
this is consistent with [2, 3]. 

The starting hypothesis of our theory is that the collective spin exci- 
tations of the doped cuprates can be described by using the proximity 
of a magnetic quantum critical point. This was proposed in Ref. [11]; 
almost simultaneously, NMR experiments in La2-jSr5Cu04 [12] showed 
crossovers which could be neatly interpreted in terms of a magnetic 
quantum critical point near a doping concentration 5 ~ 0.12 with dy- 
namic exponent z = 1. The ground state is a good superconductor at 
this value of 5, and so the magnetic transition takes place between two 
superconducting states. Evidence supporting this interpretation also ap- 
peared in neutron scattering measurements [13]. An explicit theory for 
a quantum transition between a d-wave superconductor with co-existing 
long-range spin density wave order (a SC-I-SDW state) and an ordinary 
d-wave superconductor (a SC state) was first discussed by Balents et al. 
[14] ; they focused on the case where the SDW ordering wavevector was 
exactly equal to the spacing between the two nodal points where the 
d-wave superconductor has gapless quasi-particle excitations, and stud- 
ied their role in the critical theory. However, their analysis also makes 
it clear that the nodal quasiparticles can be safely neglected for the 
generic case in which the wavevector matching condition is not satisfied 
[4], and we will mainly discuss this simpler case here. The SC-I-SDW to 
SC transition in this case is formally identical to that in an insulator, 
and the SC state has a sharp S = 1 'resonance peak' associated with 
stable S = 1 collective excitonic excitation [15]. Such a theory for the 
SC to SC-I-SDW transition was used by us [16] to predict the effects of 
Zn impurities on the resonance peak in the SC phase. 
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Figure 1. A bond-centered, coUinear SDW at the wavevector (37r/4, tt). The size 

of the arrow represents the mean spin moment, while the shading of the circle is the 
electron density. The sliding degree of freedom corresponds to a shift in the position 
of the light and dark circles. 



1. ORDER PARAMETER AND FIELD 
THEORY 

Neutron scattering experiments show that the lowest energy collective 
spin excitations at and above 5 w 0.12 reside near the wavevectors K^; = 
(37r/4, tt) and Kj, = (tt, 37r/4) (our unit of length is the square lattice 
spacing). So we write for the spin operator at the site r: 

Sa{r, r) = Re [e^^-'-^.alr, r) + e*^^--$,a(r, r)] , (1) 



where a = x,y,z extends over the directions in spin space, and ^x,ya 
are complex fields which will serve as order parameters for the SC+SDW 
to SC transition. In the SC+SDW phase, ^x,ya are condensed and the 
condensate describes the spatial modulation of the spin, by (1); in the 
SC phase, the ^x,ya are dynamically fluctuating, and its quanta are spin 
excitons. The representation (1) can describe a large variety of spin 
modulations e.g. the state with {^xa) c<. (l,i,0), {^ya) = is a spiral 
SDW in the x direction. Experimentally, however, it is clear that the 
SDW is not spiral, but coUinear, an example of a coUinear SDW has 
{^xa) Of e*^(l,0.0). {^ya) = notice that the average spin vectors on 
all sites are parallel or antiparallel. The phase 9 represents a sliding 
degree of freedom of the SDW: for the special commensurate value of 
under consideration here, the coupling to the lattice will prefer that 
9 take one of the values mr/A (a site-centered SDW) or (n -f l/2)7r/4 
(a bond-centered SDW) where n = 0...7 integer. A sketch of a bond- 
centered SDW is shown in Fig 1. Notice that the magnitude of the spin 
changes from site to site, and this implies [17] that there must be a 
corresponding modulation in the electron density, 6p{r): the latter has 
the same quantum numbers as J2a ^"ai hence we deduce that the 
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charge order can be written as 

6p{r, r) oc Re ^ [e^'^-^$L(r, r) + e^^^-^f^^lr, r)] + . . . . (2) 

a 

The period of the charge order is half that of the SDW, and the am- 
phtude of the charge order vanishes for the spiral SDW. We emphasize 
that after accounting for the long-range Coulomb interactions, the actual 
modulation in the electron charge density per site may well be unob- 
servably small. The "charge order" discussed here and in (2) should be 
interpreted in a much more general sense, as representing a modulation 
at wavevectors 2K.x,y in all spin-singlet observables which are invariant 
under time-reversal; the modulation could be larger in other observables 
like the mean kinetic energy, exchange energy, or pairing amplitude in 
the bonds between nearest-neighbor sites. 

We are interested here in the quantum transition from the SC-I-SDW 
state with {^x,ya) 7^ to the SC state with {^x,ya) = in a background 
of quiescent superconductivity. The allowed terms in the effective action 
are constrained by the underlying symmetries: these were discussed in 
some generality in [9] . Here we will be satisfied by considering a simpli- 
fied effective action which is written most easily in terms of the real and 
imaginary components of ^x,ya'- 

^xx = ^1+ ^xy = ^2 + ^xz = V3 + /g'j 

^yx = Vi + '^fW ^yy = 9^5 + if 11 ^yz = f6 + if 12 

The effective action for the real fields /j, = 1 . . . 12 is taken to be 

S^ = J d\dr{^[{drf^f + {Vrf.f + S^l] + |(^^)2| 

where a summation over the repeated fj, index is implied, and we have 
chosen units so that the velocity of spin waves is unity. Notice that the 
action 5^ has a large 0(12) symmetry of rotations in /j, space. This 
symmetry is present in all allowed terms which are quadratic in ip^, but 
is broken by a number of quartic terms [9] which are not displayed in 
(4). However, all permitted terms in the effective action do respect the 
sliding symmetry under which ^x,ya e^""^''"^^^^x,ya for integer nx,y 
The coupling s serves as the tuning parameter which measures distance 
from the quantum critical point: the SC+SDW phase will appear for 
s < Sc and the spin-singlet SC phase for s > Sc- The dynamic properties 
of the transition at s = Sc have been investigated in much detail in 
earlier work [15]. 

Now we consider the influence of the magnetic field, H, applied per- 
pendicular to the Cu02 layers. This couples most strongly to the 'back- 
ground' SC order, and so we are forced to consider the response of the 
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superconducting order parameter t/'(r). In suitable units (discussed in 
[9]), the free energy for ^(r) can be written in the famihar Ginzburg- 
Landau form 



= T / dV 



■|V'l' + ^IV'l' + l(Vr-^A)^P 



(5) 



where T is a parameter measuring the relative contributions of the mag- 
netic and superconducting energies, and Vr x A = Hz. Notice that we 
have assumed a r independent — this is permissible because the SC or- 
der is non-critical and its quantum fluctuations can be safely neglected. 

Finally, we have to couple the iJ-response of i/j{r) to the quantum 
SDW fluctuations. There arc two distinct couplings, which have rather 
different physical consequences. The first, v, is a simple coupling be- 
tween the magnitudes of the SC and SDW order parameters, chosen 
with a repulsive sign {v > 0) to account for the competition between 
these orders: 

S, = yd\dr^l{r,Tmr)\'. (6) 

Such a coupling was discussed by Zhang [18], in work focusing on the 
possibility of a first order transition between an SC phase and an insulat- 
ing phase with SDW order. The coupling v will play an important role 
in determining our phase diagram to be described below, which contains 
a second order transition between SC and SC-I-SDW phases. The second 
coupling, unlike v, recognizes the fact that the vortex lattice induced by 
H breaks translational symmetry, and so the SDW fluctuations should 
also not be invariant under the 'sliding' symmetry (the term in (6) is 
invariant under the sliding symmetry). In particular the vortex core ra- 
dius in the cuprates is only of the order of a few lattice spacings, and 
the energy of an SDW fluctuation will certainly change depending upon 
which portion of the charge order (see Fig 1) is centered on a vortex core. 
In other words, each vortex core, at a position r^,, will prefer a certain 
phase of the local charge order parameter J2a'^x,ya(VvTT). Expanding 
the charge order parameter using (3), we deduce the second term which 
couples the SC and SDW order parameters: 




The complex coupling constants Cx,y measure the pinning strength of 
the phase of the sliding charge order to some preferred value near each 
vortex core. 




Figure 2. Zero temperature phase diagram of the model Z defined in (8) and (9); 
firom Refs. [5] and [9]. The phases have the following expectation values: (i) SC: 
(v5m> = 0, (^> 7^ 0, (m) SC+SDW: {^^) 7^ 0, / 0, (m) SDW: {^^) + 0, (t/.) = 0, 
and {iv) "Normal": = 0, (V') = 0. 



We have now defined a well-posed field theoretical problem, which 
was analyzed in some detail in [9]: describe the dynamic quantum SDW 
fluctuations associated with the partition function 

Z [V'(r)] = y Vlp^,{x, t) exp {-"^ - - - S^i^ , (8) 

where the optimum value of the static SC order 'ip{r) is determined by 
the minimization of —lnZ[ip{r)] via the solution of the saddle-point 
equation 

Note the highly asymmetric treatment of the SC and SDW orders. 

2. PHASE DIAGRAM 

Our primary results for the properties of (8) and (9) are contained 
in the phase diagram as a function of s and H in Fig 2. The formal 
solution of these equations also allows solutions in which the SC order 
vanishes and ip{r) =0 everywhere: this leads to the SDW phase and the 
"Normal" phase. However we do not expect our theory to be accurate 
such a regime: quantum fluctuations of the SC order parameter are 
surely important once ip{r) becomes small, and these have been neglected 
in our theory. Our results are more precise in the small H region of 
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Fig 2, in the vicinity of the boundary between the SC+SDW and SC 
phases; indeed, the functional forms of the main results quoted below 
are expected to be exact. Upon accounting for the quantum fluctuations 
of the superconducting order, and the Berry phases associated with the 
electrons in the nearby Mott insulator [4], we expect that the "Normal" 
phase will display some sort of static charge order. 

An important prediction of our theory is the shape of the second- 
order phase boundary between the SC and SC+SDW phases at small 
H. This transition is associated with the condensation of the tp^ exci- 
ton, when the spin gap to the creation of an exciton vanishes in the SC 
phase. Detailed arguments were presented in [5, 9] showing that this con- 
densation occurs in an exciton state which is extended throughout the 
entire lattice; indeed, a variational approximation in which the exciton 
wavefunction is assumed to be simply a constant gives essentially exact 
results (as opposed to an approximation in which the exciton is strongly 
localized in the vortex cores [19]). The presence of the vortex lattice in 
ilj{r) influences the energy of this extended exciton primarily via the v 
coupling in S^. When spatially averaged over r, the dominant change 
in the average value of |V'(r)p, and hence in the energy of the exciton, 
arises from slight suppression of superconductivity in the superflow re- 
gion surrounding each vortex core. The much smaller vortex core region 
always has a significantly weaker effect on the exciton energy. It is a 
simple matter to compute the correction to the exciton energy from the 
superflow. The average kinetic energy of the superflow is known from 
standard Ginzburg-Landau theory to be ~ (H/Hc2)lii-{Hc2/H) (here 
Hc2 is the upper critical field at which superconductivity disappears), 
and the coupling v in (6) therefore leads to a corresponding change in 
the efi^ective value of s controlling the exciton energy: 



where C is a positive constant. The critical flcld at which the spin gap 
vanishes is therefore determined by s^g = Sc, and this leads to our result 
[5, 9] for the phase boundary between the SC and SC-I-SDW phases: 



Note that the phase boundary approaches the H = limit with vanishing 
slope: consequently a relatively small field applied to the superconductor 
for s > Sc will drive the system into the SC+SDW phase. This is our 
explanation for the shift in the energy of the dynamic spin fluctuations 
seen in [20]. 




(10) 
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The phase diagram of Fig 2 leads to a number of predictions [5, 6, 9] 
for observables in the SC and SC+SDW phases, some of which have been 
tested in recent experiments [1, 2, 3]. We discuss theory and experiment 
in the two phases in turn in the fohowing subsections. 



2.1. STATIC CHARGE ORDER IN THE SC 
PHASE 



The SC phase has {ip^) = 0, and so the SDW fluctuations are dy- 
namic and the spin exciton only exists above a finite energy gap A. 
The superconducting order V'(r) is suppressed in the vortex cores, and 
so this region should exhibit characteristics of the doped spin-gap {i.e. 
paramagnetic) Mott insulator, as was argued in [21]. Paramagnetic 
Mott insulators, and their response to doping with mobile charge car- 
riers, were studied at some length in [4]: it was argued that the most 
likely candidates had bond-centered charge order which survived in a 
superconducting state for a finite range of doping — this work will be 
reviewed further in Section 3. Reasoning in this manner, Ref. [6] pre- 
dicted that static charge order should appear in and around the vortex 
cores, coexisting with dynamic spin fluctuations in the SC state. Order 
of this type appears to have been seen in the recent STM experiment on 
Bi2Sr2CaCu208+5 [1]. Our approach should be contrasted from other 
recent investigations [22] which have both static spin and charge order 
in the vortex core. 

The spatial extent of the charge order was computed in the field theory 
in (8), (9) in [9]. The gapped spin exciton, (p^, views the vortex lattice as 
a periodic potential, and consequently its dispersion develops the Bloch 
structure of a particle moving in a periodic potential. To zeroth order 
in the pinning terms, C,x,yi the two-point 4>fj, Green's function is diagonal 
in the /x index, and its diagonal component can be written as [5, 9] 



where k is a Bloch momentum which extends over the first Brillouin 
zone of the vortex lattice, /x is a 'band' index, S^k(r) are the Bloch 
states (H^k(r-|-R„) = e*'^'^"H^k(r) where R^, is an vector connecting two 
vortex centers), -E^(k) are the energy dispersions of the various bands, 
and Un is an imaginary Matsubara frequency. The lowest energy exciton 
has energy -Eo(O) = A and wavefunction E!oo(r); this wavefunction is 
sketched for typical parameter values in Fig 3. 

Now let us consider the influence of the pinning term in (7). Combin- 
ing (7) with (2) it is simple to see that to first order in the C,x,y, static 




(12) 
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Figure 3. Spatial structure of the spin and charge correlations in the SC phase at 
two different magnetic fields (squares and triangles; the squares have the lower field). 
The vortices are centered at a; = 0, f , 2. The orientation of the spins fiuctuates with a 
frequency of order the inverse spin gap h,/ IS.. The open symbols represent the lowest 
energy exciton wavefunction Hoo(r). The nature of the microscopic spin correlations 
can be understood by recalling that this wavefunction is the envelope of the order 
Fig 1. After including the pinning term <Spin, static charge order is induced, and 
its envelope f2(r) (defined in (13,14)) is shown above. The numerical results were 
obtained by Ying Zhang and reported in [9]. 
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charge order appears in the SC phase, with 



{Sp{r)) cx Re C^e' 



,i2Ka;-r 



(13) 



where 



(14) 



A plot of the function fl{r) is sketched in Fig 3, along with the corre- 
sponding Hoo(r)- At the higher field (triangles), the spin exciton wave- 
function Soo(r) is essentially constant across the entire system, with only 
a weak modulation induced by the vortex lattice; nevertheless, at the 
same field the charge order fi(r) has a strong modulation on the scale of 
c/(2A), where c is a spin- wave velocity [8]. At the lower field (squares), 
there is larger modulation in the spin exciton (but Hoo(r) only decays to 
half its maximum value) , and again the decay length of the charge corre- 
lations is about half that of the spin correlations. The spatial form of the 
lower field ^{r) in Fig 3 is quite similar to envelope of the modulation 
observed in [1]. 

The STM experiments of [1] actually measure the modulation in the 
local electronic density of states (LDOS) in a range of energies as a 
function position in vortex lattice. A great deal of information is, in 
principle, contained in the energy and spatial dependence of the LDOS 
modulations. We have recently analyzed [8, 10] simple models for the 
coupling of the electronic quasiparticles to the collective SDW and SC 
degrees of freedom that have been discussed here: these lead to predic- 
tions for the LDOS modulations, which can be usefully compared with 
the STM data — the reader is referred to the papers for details. 

We also mention here the recent STM experiments of Howald et al. 
[23] which have observed charge order similar to that of [1] but in zero 
applied magnetic field; this order has (presumably) been pinned by im- 
purities. 



The SC-FSDW phase has (c/?^) ^ 0, and hence via (1), (2), and (3), 
there is both static spin and charge order. Neutron scattering measure- 
ments have so far only succeeded in observing the static spin order, and 
so we will restrict our discussion here to the spin moment. 

In the presence of an applied magnetic field, the vortex lattice will 
spatially modulate the value of (0^(r)), and this should, in principle, 
lead to satellite elastic peaks [19, 5, 9] surrounding the main elastic 
Bragg peaks at iK^,, ±Kj, observed in neutron scattering. However, 



2.2. 



STATIC SPIN MOMENT IN THE 
SC+SDW PHASE 
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Figure 4- Spatial form of {ipfj,{r) in the SC+SDW phase. Notice the vertical scale, 
which shows that the overall modulation in the size of the order parameter is quite 
small. The numerical results were obtained by Ying Zhang and reported in [9]. 

our discussion above on the dominance of the superflow effect shows that 
this modulation occurs predominantly on the scale of the vortex lattice 
spacing [5, 9], and not on the scale of the vortex core. The resulting 
satellite peaks are consequently found to be extremely weak, as will be 
clear from our results below. 

We show a sketch of the spatial form of (92^ (r) at a point in the 
SC+SDW phase in Fig 4. All the coupling constants in the theory are 
the same as those used for the results in Fig 3 for the SC phase. Only 
the parameters s and H are tuned to move the system between the SC 
and SC+SDW phases. The spatial Fourier transform of Fig 4 determines 
the strength of the elastic Bragg peaks that will be observed in neutron 
scattering. In particular, the dynamic structure factor of y?^ has the 
form 

5^(k,a;) = (2vr)5H^|/Gp(27r)25(k-G) (15) 

G 

where G are the reciprocal lattice vectors of the vortex lattice. The 
reader should keep in mind, via (1) and (3), that the experimental 
structure factor is obtained from (15) by measuring wavevectors from 
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Figure 5. Magnitudes of the elastic scattering peaks in (15) obtained by the spatial 
Fourier transform of results like those in Fig 4. The field H is measured in units of 
H^2j the value of the critical field at which all four phases meet in Fig 2. The plot 
(a) is at s = Sc, while (b) is for s < Sc- All other couplings are identical to those in 
Figs 3 and 4. The numerical results were obtained by Ying Zhang and reported in 
[9]. 

the SDW ordering wavevectors ±K^, ±Ky. We show typical results for 
the field dependence of the strength of the central peak, |/oP, and also 
for the first satellite peak, |/giP (where Gi is the smallest non-zero re- 
ciprocal lattice vector) in Fig 5. The strong observable effect is in the H 
dependence of the central peak |/oP- The same superflow effects which 
were responsible for (11), also dominate in determining the average mag- 
nitude of the SDW order in the SC+SDW phase: we can estimate the 
enhancement of magnetic order by the superflow by assuming that /o is 
determined by Se//; and then (10) leads to 

\fo\\H)-\fo\\0)^Hlnil/H). (16) 

the lines in Fig 5 are fits of the full numerical solution to (16). Arovas et 
al. [19] had discussed nucleation of static magnetic order in the vortex 
core in what was an SC phase. In contrast, we claim that there are no 
static spins in the vortices in the SC phase, and only pinned static charge 
order around each vortex as discussed in Section 2.1. The static moments 
appear only when there is bulk magnetic order as in the SC+SDW phase 
(as pointed out in [5, 9]), and here it is the contribution of the superflow 
which always dominates leading to (16). The linear small H dependence 
of the ordered moment proposed by Arovas et al. is not valid in either 
the SC or the SC+SDW phases. 
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The results in Fig 5 also compare well with the experimental obser- 
vations in the overall scale of both the field and the magnetic moment: 
those in (a) are quite similar to the results of [2], while (b) matches well 
with [3]. 

Finally, note that Fig 5 shows that the satellite peaks are unobserv- 
ably small. This is related to relatively slow modulation induced by the 
superflow in the moment in Fig 4. Ref. [9] predicted that the influence of 
the vortex lattice may be more easily observable in the dynamic exciton 
band structure in the SC phase. 

3. WHY DOES THE CHARGE ORDER HAVE 
PERIOD 4 ? 

The SDW (Kx,y) and charge {2'Kx,y) ordering wavevectors have so 
far been arbitrary parameters in our phenomenological theory, and the 
structure of this theory is largely independent of the values of 'Kx^y (some 
high order terms in the action are permitted only for certain commen- 
surate values of i^x,y)- The determination of ^x,y requires, instead, a 
lattice scale theory of the doped antiferromagnet. 

The mechanism of the charge-ordering instability in doped antifer- 
romagnets has been discussed using a number of different theoretical 
perspectives by several other workers [24, 25, 26, 27, 28]; in general, 
many possible charge-ordering periods emerge in their works, and there 
appear to be no fundamental principles restricting the possible periods, 
or whether the charge-ordering is site or bond centered. Here, we will 
briefly recall our theoretical work on charge ordering [4]: its point of 
departure is the theory of the magnetic quantum critical point in the 
Mott insulator. Our perspective led to bond-centered charge ordered 
states, without long-range magnetic order, and with co-existing d-wave- 
like superconductivity. The predicted evolution of the wavevector of 
this ordering, as a function of hole concentration, is shown in Fig 6. 
Note that the period, p, is always pinned to be an even number i.e. 
2K-r = (27r/a)(l/p, 0), where a is the square lattice spacing. There is a 
large range of 5 values in Fig 6 where the period is pinned at p = 4: this 
corresponds to the values of 2'Kx,y observed in the STM experiments 
[1, 23]. Both experiments also see the modulation appearing simultane- 
ously at 2K.X and at 2K.y, leading to a checkerboard appearance in real 
space. Such checkerboard patterns have been considered previously for 
the p = 2 case [29, 30, 31, 32], and were found to have an energy very 
close to that of the state with charge order only along a single direction; 
we can expect that a similar result applies for p = 4. 
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Figure 6. Evolution of the charge ordering wavevector upon doping a paramagnetic 
Mott insulator, from Ref. [4]. The wavevector is measured in units of 2n/a. At 5 = 0, 
the Mott insulator has period 2, which is associated with the appearance of bond 
order (see Fig 7). Beyond a range of very small 5 values, the ground state is also a 
superconductor. Full square lattice symmetry is restored above 6 « 0.175, when the 
ground state becomes an ordinary d-wave superconductor. 



Further support for Fig 6 has emerged from recent neutron scattering 
measurements of Mook et al. [33] in the under-doped superconductor 
YBa2Cu306.35: they observed static charge order, and dynamic spin 
correlations, with the the charge order period pinned rather precisely at 
j» = 8. This is in accord with the plateau at p = 8 over a range of small 5 
values in Fig 6, and should be contrasted with the continuous evolution 
of ordering wavevector with 5 which is usually assumed in the "Yamada 
plot" [34, 35] for periods larger than p = A. 

We conclude this paper by briefly recalling the physical ingredients 
behind the results in Fig 6. A review of these arguments has already 
been presented in [7] , and we present here a synopsis in Fig 7. In moving 
from an undoped Mott insulator, like La2Cu04, to a high temperature 
superconductor at optimal doping, at least two quantum phase tran- 
sitions must take place: one involving the loss of magnetic order, and 
the other the onset of superconductivity. Theoretically, it is very useful 
to disentangle the two transitions by imagining that we have a second 
tuning parameter at our disposal, in addition to the doping 5. We use 
this second parameter to first destroy the magnetic order in La2Cu04 
while remaining at 5 = — a specific possibility for such a parameter 
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Further neighbor exchange 
with square lattice symmetry 




Figure 7. Schematic T — phase diagram (from Refs. [4, 7]) for the high tem- 
perature superconductors as a function of a ratio of the near neighbor exchange in- 
teractions and the hole concentration, S; e.g. the vertical axis could be J2/J1, the 
ratio of the first to second neighbor exchange. The shaded region has charge order. 
The hatched region has broken spin-rotation symmetry with (5) 7^ 0, and at least 
one of the order parameters ^x,ya is non-zero; (S) = elsewhere. The unfrustrated, 
insulating antiferromagnet with long-range Neel order is indicated by the filled circle. 
At 5 = 0, there is an onset of charge order above the point A, while spin-rotation 
invariance is restored above B. The nature of the charge orders as determined by the 
computations of Ref. [4] are indicated at the top of the figure; numerous transitions, 
within the gray shaded region, in the nature of the charge ordering are not shown. 
The ground state at very low non-zero doping is an insulating Wigner crystal and 
there is subsequently a insulator-to-superconductor transition; superconductivity is 
present over the bulk of the 5 > region. The central idea behind our approach is 
that many essential aspects the spin excitation spectrum of the insulating, paramag- 
netic region {5 = 0, (5*) — 0) lead to a simple and natural description of the analogous 
properties of the d-wave superconductor. 
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is a frustrating second neighbor exchange interaction (see Fig 7). De- 
tailed arguments have been given (for a recent review see [36]) that the 
paramagnetic Mott insulator so obtained has bond-centered charge or- 
der, as indicated in Figs 6 and 7. The second theoretical step of doping 
the paramagnetic Mott insulator can be reliably addressed in a large N 
theory [4], and one eventually obtains a d-wave superconductor which 
fully respects the symmetries of spin rotations and lattice translations 
and rotations. In this approach, the charge order of the superconductor 
without magnetic long range order evolves from that in the paramag- 
netic Mott insulator. Further, the Cooper pairing in the superconductor 
is also connected to the singlet electron pairing present in the bond- 
ordered paramagnetic Mott insulator. These connections led us to the 
results in Fig 6, and to the proposal of charge order nucleation by vor- 
tices in [6]. 
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